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I. INTRODUCTION 


Many important issues in the control of Urge space structures are intimately 
reUted to the fundamental problem of parameter identification. Very often, a 
complicated structure can be adequately modeled for certain operations by the 
fitting of a rather simple model with a number of free parameters. This simple 
model then can be referenced for necessary control operations. Important 
applications include the many space station designs which are based on the assembly 
and Joining of discrete module* by crew memoers. This crew-assisted construction 
will result in a configuration * oich is a Urge-scale composite of many structural 
elements and whose static and dynamic characteristics cannot be adequately 
modeled in advance. In fact, any modeling will require periodic updating as more 
modules are added to the system and as the structural properties of the elements 
slowly change over the lifetime of the station. 

One might also ask how well this identification process can be carried out in the 
presence of noisy data since no sensor system U perfect. With these considerations 
in mind our algorithms are designed to treat both the case of uncertainties in the 
modeling and uncertainties in the data. 

This paper serves as a companion to [6] where the analytical aspects of 
maximum likelihood identification are considered in some detail. Here we focus on 
the questions relevant to the implementation of these schemes, particuUrly as they 
apply to models of Urge space structures. Our emphatU will be *n the influence of 
the infinite -dimensional character of the problem on finite- dimensional 
implementations of the algorithms. We highlight those areas of currer . and future 
analysis which indicate the intezpUy between error analysU and possibU truncations 
of the state and parameter spaces. 

2. MODELS 


As in (6], we consider the systems of the form 
A(©)u(e) » a B(6)<i> + C(©)f 

<j) 

( 2 . 1 ) 

y(0) = H(6)u(6) to^T) 

Here A is a formally self-adjoint elliptic differentUl operator defined over the 
spatUl domain ft; the integral operator <t> is reUted to A by 

A<t> - I , (2.2) 
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where I is the identity. B end C are appropriately dimensional operators that model 
the influence of the process error o and the input f on the state u. H is an 
operator that characterizes the state-to-observation mapi <■> amd n are model 
errors that form the model error vector 

e = (o, T|]t (2. 3) 

and f is a deterministic input. Conceptually the error vector £ represents spatial 
white noise and is characterized by the covariance operator 

E(££*) - I. 

a and cr are non- negative scalar weighting parameters that respectively 

u) T| 

measure the relative importance of the modeling error and the measurement error. 
Thus, the limit cr -* 0 corresponds to the case of perfect modeling while the limit 

-* 0 corresponds to the case of perfect measurements. 

6 is the possibly infinite-dimensional parameter which must be estimated. For 
simplicity we shall generally consider cases where the parametez dependence is 
restricted to the operator A. Furthermore, we assume as in [6] that the parameter 
enters linearly into the expression for the potential energy of the system. Thus we 
assume 


A(0hi * D*(6Du) 


(2.4) 


where D* denotes the formal adjoint of Dj the corresponding potential energy is 
given in terms of the appropriate state-space inner-product: 

<A(6)u,u> » <8Du,Du>. (2.~, 

And finally, the deterministic and stochastic forcings will be localized to 
discrete points which might correspond to actuator locatir's. Similarly, the 
observation map returns a vector of observations at discrete points which might 
correspond to sensor locations. We assumed that there are N point-sensors at 

locations {f ^ and point- actuations at locations {f 


Because of these last assumptions, many of the relevant calculations outlined in 
[6] reduce to matrix and vector manipulations, hi this paper the notation X will 

refer to a finite-dimensional vector whose k-th component is given by g^\ 

Similarly G is the notation for a matrix whose (i,))- component .is given by The 

relevant dimensions of vector and matrix quantities will always be clear from the 
context. 


After taking formal limits in the system (2.1) we have: 

Au * a B^ ♦ cT 

a 

f - Hu + zf 


( 2 . 6 ) 
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where 


- f AN >VT 

o * (tt ,...» o a } 

(2.7a) 

- - Ci) (NAT 

T) s (T1 5 ) 

(2.7b) 

Hu(x) ■ (u(£ u(fj^ ) T 

(2.7c) 

s 


N , 


H ♦T*E6(x-C i )g (i) 

(2.7d) 

i=l 


N 

• 


cF= £ 6(1-^) f (l) 


i=l 

(2.7e) 

C*u(x) » (u(f i ) u(J N )) T 

(2.7f) 

* 


N . 


Bo = l 6 (x-^) o (l) 

(2.7g) 

i=l 


B*u(x) - (u(f y u(f N )) T 

(2.7h) 


a 


*m in the cate of more general sensing and actuating systems, the modeling 
requirements for the system can be reduced to solving equations of the form 

Au * f. (2.8) 


Thus the discussion in this section will focus on how the infinite-dimensional 
structure of the system (2.8) influences the choice of finite-dimensional 
approximations which can be made. In this paper we consider two specific 
structural models: a string under tension and a wrap-rib antenru. 


Let £ be a distance coordinate measured in meters along a string of Length L which 
is also given in meters. Let u(x) be the displacement in meters and let <£(x) the 
tension parameter given in units of newtons. The forcing density is given by f(x) in 
units of newtons/meter. Then the energy potential [3> is yiven by 


A A 

V(u) 


if 


-1 


a(xXu'(x)) a d£ - f(r) u(x) d£. 


(2.9) 
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Of course the energy potential is given in units of newton- meters. The equations of 
motion can be derived immediately based on the principles of the calculus of 
variations but it will be convenient to first transform to dimensionless coordinates. 
Let a, be some characteristic value of the tension parameter. We introduce the 

dimensionless variables: 

A 

X 

X = 

L 

a<x) 

a(x) . 


u(x) = 


A A 

u(x) 


L 


aa 


f(x)L 


f(x) = — 

A 

a* 


( 2 . 10 ) 


A A 

V(u) 


V(u) = — 

A 

a*L 

and the potential expression becomes 


V(u) 


•j" f a(x) (u'(x)) a dx - f 

Jo Jo 


f(x)u(x)dx. 


( 2 . 11 ) 


Por simplicity we prescribe boundary conditions corresponding to fixed end points: 

u(o) >u(l)>o (2.12) 

Then arguments based on the calculus of variations give the system 
(a(x)u'(x))' - f(x) 

( 2 . 13 ) 

u(o) » u(l) * 0 0 < x < 1 

This example has been studied many times in the classical literature but an 
analogous approach gives comparable expressions for much more complex systems. 

We consider now a planar model for a wrap-rib antenna which is used to study 
out-of-plane vibrations (see Figure 1). The antenna model comprises N gores 
(subsections) modeled by ‘nterconnected ribs and mesh. Since the transformations 
are similar to those used in the case of the string, we immediately write the 
potential expression with dimensionless coordinates. Let the vector of rib 

displacements be u(r) where the k-th component of vTis the displacement of the 
k-th rib (0 < r < 1). Let the vector of mesh displacements be"?tr,6) where the k-th 
component of Vis the displacement of the k-th mesh sector (0 < r < 1, 0 < 6 < 1). 
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Fig. 1 Simplified model for wrap-rib antenna 

Based on analysis of actual antenna designs, our model consists of N identical 
beams fixed at a central bub. Stretched between the beams are N identical 
anisotropic membranes. The potential equation is given by 



f 1 - * _ fl f 1 - * _ 

- I F D u dr - I I F„ vrdrd© 

Jo R J 0 Jo " 

Here the coefficients {G^} are related to the physical parameters of the beams 
and membranes thusly: 


El 

o 

V — 

oL 

T e L a 
r o 

G = (2.15) 

a 

T L a 

T e L 

V 

o9 

o 
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E and I are respectively the Young's modulus and the moment of inertia of the 
beams. T and Tq are respectively the radial and circumferential tensions of the 
membrane. L is the radius of the antenna and 9 q is the angular width of a sector; 
that is, we have 


2w 

8 = , ( 2 . 16 ) 

o 

N 

where N is the number of gores. Finally, a is some convenient scaling parameter 
with the dimensions of energy (nt-m). We note that the physical forcing densities 

and Fjg having respective dimensions nt/m and nt/m* were rescaled according to 


- L 7T 

^R = _ P R 


F M = “ P M 


( 2 . 17 ) 


Appropriate geometrical boundary conditions follow from fixing the center and 
attaching each of the ribs to its adjoining membranes: 


u| 


dr 


u| 


r=o 


r*o 


( 2 . 18 ) 


7 1 = C v| j u 

9=0 9=1 

Here C is an N x N periodic matrix: 

'0 1 


( 2 . 19 ) 


. 1 
0 
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As in the case of the string, the equations now follow from arguments based on 
the calculus of variations: 




(2.20a) 


1 


r 



(2.20b) 


with the additional natural boi*ndary conditions: 


8*\T 

dr* 


« 0 


r=l 


3*u 

** r=l 


dv* 

”5T 



( 2 . 21 ) 


One of the focuses of this study is the consideration of nonmodal approaches in 
the finite -dimensional approximation schemes. In practice, this generally will mean 
directly solving a linear system of equations rather than proceeding from some 
finite modal synthesis. But the infinite -dimensional structure of the system (2.8) 
also can influence the particular finite-dimensional approximation schemes used. 
Our approach is sufficiently general so that any adequate finite e 1 lent model of 
the system (2.8) should yield adequate numerical approximations. Bv one can often 
do much better for a particular model or a particular class of models. 


We use the antenna model to illustrate the point and make some observations 
that should influence the approximation schemes regardless of which finite element 
or finite difference scheme is employed. We emphasize that these considerations 
also apply to much more complicated antenna models which share salient features 
with the system (2.18) - (2.21). First we note that the structure is periodic in the 
6-direction. This cyclic symmetry leads to considerable savings in the computation 
of solutions to (2.8). This can be deduced from either the differential equation or 
the energy expression (2.14). The periodic matrix C can be diagonalized by means of 
a finite Fourier transform [1]. That is, let U be the N x N matrix whose (j,k) 
component has the form ~ 


0,k) 

U 


exp (i -^j — 0-1) (k-1)) 
J N 


( 2 . 22 ) 
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We then have: 


U*U * I 

<">w esv '"W» 

U*CU = A = diag exp (i (j-1)) (2.23) 

~ ~~ ~ N 

j 

This transformation decouples the system since the potential expression (2.10) 

with 


U*u 


7 = U*v 

IT = U*F 
*R ~ r R 

F., = U*F 

M ~ m 


(2.24) 


has the same form as the original system except that the matrix C is replaced by 
the diagnonal matrix A. The differential system (2.20) is likewise transformed. 
Thus any particular solution of (2.8) can be expressed in terms of N subsystems each 
comprising a single rib counled to a single membrane. Since the cost of solving a 

m - dimensions 1 linear system is 0(m*) this represents a considerable computational 
savings. 


The balancing of terms in the equation also can influence the choice of 
discretization. Based on a report by Lockheed on the specifications for a 55-meter 
wrap-rib antenna with 48 ribs [2], the following nominal parameter ranges were 
derive : 


L ~ 27.5m 

6 ~ 1.31 1(T X 
o 

I ~ 1.31 1(T* m 4 

(2-25) 

B ~ 9.72 10 l ° nt/m* 

T r ~ 1.75 10~ x nt/m 

T 0 ~ 3.50 10' 1 nt/m 

This lives the proper scalings in system (2.20). For simplicity we take 

a - T, L */Q and we have 
G o 
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2.29 


G ~ 

1 

G ~ 8.58 10“’ (2.26) 

2 

G ~ 1 

s 

This nesni that the radial terms of the mesh potential are comparatively small 
except when the radial derivatives are large. How this affects the structure of the 
system is demonstrated in the following example (see Figure 2). 



Fig. 2 Single antenna gore 


Example: 


1 3 / 3v\ 1 

T ~5F \~SFf ~T 


_3fv 

30 a 


v (r,0) 

0*o 

»f (r) 
i 

o < r < 1 
o < 0 < 1 

(2.27) 

v (r,0) 

0*1 

* f (r) 
2 

(f (o)*f (o) * o) 
1 2 


5T (r,6) | 

* 0 




r»l 


hi example we study the equations for a sector of membrane where the 
prescribed boundary conditions depend on the adjoining rib displacements (Mr) and 
Mr)). For simplicity we take the forcing on the mesh to be aero although the more 
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general case coaid be handled in a similar fashion. We are of course interested in 
the case where 


0 < e « 1 


(2.28) 


which corresponds to the parameter ranges (2.26) in (2.20). Physically one expects 
that the radial terms contribute little to the static behavior except perhaps at the 
boundary where the gradients may become large. One is also interested in the 
behavior near the comers (r,0) * (1,0) and (r,6) * (1,1) since some «faignUr behavior 
may be possible. Using the techniques of singular perturbations, (see for example 
[4]), one can show that as e approaches zero we have 

v(r,0) = f (r) (1-0) ♦ f (r)0 
1 2 


♦ e V * Z ( - f Vl> ♦ f ' (l>(-lf ) 

♦ 0(e) 


2 

(nrr ) 2 


exp 



sin (mrO) 

(2.29) 


This expansion could be continued to higher orders, and, as noted before, a more 
complicated expression would result from forcings on the membrane. One possible 
approach to the numerical solution of the system (2.20) would by the elimination of 
the mesh behavior entirely be substituting an expression similar to (2.29) into the 
beam equations (2.20a). Then one would have only equations along the beams to 
solve. If higher-order accuracy on the mesh is required, one could then apply 
finite-element techniques to the system obtained after linearization about the 
asymptotic expansion for the mesh behavior. Finally we note the appearance of 
logarithmic singularities in the mesh gradients (dv/dv, dv/90) as the comers (r,0) = 
(1,0) and (r,9) * (1,1) are approached from the interior of the mesh. This 
consideration should also influence any finite element approximation of the mesh. 

We emphasize that this analysis applies not only to the simplified antenna 
model we have considered but would hold for more elaborate configurations where a 
similar structural balance of terms governs the system. Thus, many modeling 
options can be considered for parameter identification in important classes of 
structures if one does not insist on a traditional modal characterization of the 
system. 

3. THE LIKELIHOOD FUNCTIONAL 

A detailed discussion of the likelihood principle is given in [6]. The functional 
we consider is the negative logarithm of the likelihood ratio associated with the 
detection of a Gaussian signal in additive Gaussian noise; this framework is 
traditional in the theory of communication and signal detection. 

In accordance with the discussion given in [6], and the notation discussed in 
Section 2 ,the log-likelihood functional is given by, 
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(3.1) 


J(6,y> Tr (log loJjL + o* R]) 

♦-^(y-n5(0))* [a 2 l* R] “‘Cy 1 ®©)) 

m TT^ 

“a (1A V*^ 


where the expected mean and the covariance operator are given by: 


m(0) * H(0) $>(©) C(0) f 

(3.2) 

R(0) = H(0)<i>(e)B(e)B*(e)d>*(0)H*(e) 

Here y^*y^ indicates the Buclidian inner product in the N^-dimenaional apace to 
which the observators belong. 

From the assumptions of Section 2, it is easy to see that R is an N^x N ^ matrix 
whose (i,j) component is given by 


R (U) 


k=l 


(3.3) 


N 

k-l 


where g(x|£) is the point-source solution of the underlying elliptic system 

Au * 6(x - {) (3.4) 

with the appropriate boundary conditions. And likewise the expected observation 
has the form 

55 ,= G f 

r-(f (l> W) T (3.5) 

G (i,,) a ^CjICj) 

where G is an N xN matrix. 

~ 8 A 

We note that (3.1) also differs from (1.3) of [6] in accordance with the 
introduction of the positive weighting parameters, a ^ and a^, Into the system (2.1). 
An equivalent form for the likelihood functional follows from a rearrangement of 
terms. 
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(3.6) 


J(0,yj = *2 Tr log [I^+ uR] 

♦ (1 /a^ a (y^3(0))* [I^+ vR] -1 (y^m(0) 
- (1/cM 2 y*y 


where 


V - 



(3.7) 


This form is useful since one can arrive at the functional given in [6] directly by 
the substitutions 

A 

R ■* ii R = R 

Aw 

y -» (1 ~ y (3.8) 

in -* U/a ) m * m 
ti 


This correspondence allows one to use the algorithms derived in [6] directly on the 
functional 


J (0,y) =4 Tr Log l]U R] 

+4 (y - *b( 0))* [i + £] -1 (y - m) 


1 a A 

- T y*y 


(3.9) 


The goal is to find the parameter value 0 which minimizes the log-likelihood 
functional; that is, we wish to solve 

min 1(0, y) , (3.10) 

0 


where 0 ranges over some appropriate infinite-dimensional space. Assuming that 
the functional hr.s a Prechet derivative and satisfies an appropriate convexity 
condition, one can restate the problem (3.10) as 

81/30 (0,y) « 0 (3.11) 

Both problems (3.10) and (3.11) have been studied in a variety of contexts (See, for 
example, [5]). 

Since both the parameter space and the state space are infinite-dimensional, 
one must make dual approximations in order to achieve problems that are 
finite -dimensional and therefore computationally tractable. Thus in practice one 
solves a sequence of problems of the form 
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(3.12) 


AAA 

min J(0,y) 


£ 


or 

af/ae (£,y) - o 


(3.13) 


where the state-space and parameter- space have been replaced by finite 
dimensional spaces. Then the problem reduces to a finite minimi sation problem 
which can be treated numerically by a variety of techniques (see, for example, [1,5]). 


The state-space can be approximated by a finite -element space which is 
appropriate for approximating solutions to (3.4), and the parameter-space can be 
conveniently represented by a spline-based space. Let be the dimension of the 

finite-dimensional approximation to the state-space and let Nq be the dimensions of 
the finite-dimensional parameter space. This leads to the natural substitutions 


N 


u -» £ u 
k=l 


(k) 


♦k ( *> 


N e 

e -» i e (k> <c k (6) 
k»l 


0 a [6 


,(i) 




(3.14) 


U a [U 


(1) 




> 


where the sets and {r^} give the basis elements for the state and parameter 
spaces respectively. 

It will also be convenient to consider the state-space inner-product with a 
weighting given by the basis elements of the parameter space. Thus we define 


< U,V = <U, lfj(0)V> 
0 • {!»...» Nq}) 


(3.15) 


In the following we restrict our attention to these finite-dimensional problems, 
and, when the context is clear, we suppress the ''-notation. Questions concerning 
the convergence of the numerical schemes and the general relationship between the 
infinite -dimensional and finite -dimensional problems will be discussed more fully in 
a future report. 
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Most nonlinear optimization techniques require solving linearized systems 
iteratively, and consequently one must solve systems of the form (3.4), where the 
dimension may be quite large. Since the complexity of solving an m-dimensional 

linear system is 0(m*), the speed of convergence of the iterates is an important 
consideration. With this in mind, we emphasize the use of quasi-Newton methods 
for the solution of (3.13). Consequently much of the resulting effort is directed 
towards deriving adequate approximations for the Nq- dimensional Jacobian vector 

3J/30 and the Nq x Nq Hessian matrix 3 a J/30 2 . 

We briefly outline the procedure here; as noted previously, a more complete 
description is given in [6]. In general for the finite problems, the dimension of the 
state space (N^) is much larger than the dimension of the parameter space (Nq), the 

number of sensors (N ) or the number of actuators (N ), and so it is preferable to 

S A 

carry out the necessary manipulations in spaces whose dimensions do not depend on 
the dimension of the state space. 


Therefore, as in [6] we represent calculations in terms of the eigen-structure of 

the N x N matrix R. 
s s ~ 



= X 



X k ‘ M a k <° 5 a r < 7 ) 


(3.16) 


From the spectral components of R we define useful quantities as given in [6]. 

From (3.2) we have the expected observation 

m = Hd>c7 

(3.17) 

“k * V 5 


and also we define the filtered observation 


"z" = Ly* + (I - L) m 

~ ~ ~ (3.18) 

\ ■ V r 

where the N x N matrix L is given by 

S3 ^ 

_y. 

L = I - (I + R) 

~ ~ ~ ~ (3.19) 

= T. (1 - COSG^) 

k 

and the related matrix K is given by 
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(3.20) 


K - a ♦ R) Vl - I 

r>J fKj ^ 

• £ (seca^ - 1) <? k ^ k * 
k 


For algebraic convenience we also define the residual of the process: 
"e = y-7 

e k 3 V r 

The gradient of L is represented by 

3L/ae (,> - z E.’ n AV 

k m#k 


where the coefficients {aj^} are given by: 


- (sin a^)* <Dp k> Dj^ >. 


, k « m 


‘km 


((X^ X k ) / (\£ - \^)) (cosc^ - cosa m ) • 

[x k <D »m' ^ V X m <Dp k' > )> 


, k 4 m 

J 


For later convenience we derive another form for the coefficients {a^ |y }. 
standard trigonometric identities one can easily verify the relation 


cosa. - cosa (cosa .) a (cosa ) a 

* m k m 

(tana. ) a - (tana ) a cosa + cosa. 

km tn k 

This leads to an alternate form for the coefficients 

- ( coso^* t(X k >* <Dp fc , > } ] , k-m 

- (cosa k cosa m > a / (coso^ + cosa m > . 

IX m V J <Dp m’ ^ V H K m <D *V ^m >f ’ k * m 



(3.21) 


(3.22) 


(3.23) 


Using 


(3.24) 


(3.2S) 
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The point of this last derivation is that the bracketed terms reduce to simpler 
expressions. From (2.2), (2.7), and (3.4) one can easily show: 


*k . X-‘ «* H* * k 
N 


= x ’ 1 1 g (xie^y 

and also from (2.7) we have: 


0> 


x = B*p 
m m m 


N N 

a s 


(j) 


k.l i-t . 


By this we have: 

^m^k'^mV^k^m 

where the N x N matrix A. has the form 
s s ~) 

N a 

Aj k,m) *Z g tf^) <Dg (x|f k ), Dg (x|C t ) >j • 

1>1 

And similarly we have the useful relation 
H <DP k . DEyVii 7 
u * 4>Cf 

where the N xN matrix B. is given by 

8 A 

B (k,m) = <Dg (x(Ck)i Dg (x,f m »j 

We now give expressions for the gradient and the Hessian in terms 
quantities given above. As in [6] the gradient can be represented as 

3J/ae (i) » -£ sin* tan <Dr k . 

-l t y (dz v /30 (i) ) 

Here the spectral coefficients 3z k /3e (i) are given by 


(3.26) 


(3.27) 

(3.28) 

(3.29) 

(3.3b> 

(3.31) 
of the 

(3.32) 
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(3.33) 


a V ae< ” * 1 CMa m * L e k 

m 

- (cosa^) CDpj,, DS>jJ 

Exact expressions for the Hessian are given in [6]s in general, however, all 
terms need not be estimated to give an adequate approximation. In particular, the 
calculations are much simpler if the terms with second-order derivatives can be 
ignored. The simplest approximation comes from only keeping those terms which 
contribute to the expected vt of the Hessian. Thus, from (1.12) of [6] the 
(i,j) - component of the N x N Hessian approximation M is given by 

8 S < “ w 


* tr i3L/ae (1) (i + K) 3L/ae (J) a ♦ K) i 


(3z?30 (i) ) * (3z736 (j) ) 

(3.34) 

* tr tXi Xj] ♦ E eyae (i) dz k /3e (}) 


where by (3.20) and (3.22) we have vj m,k> - aj^seco^. 



This estimate is justified when the covariance is small as one might expect if the 
number of measurements is large. This point will be investigated more rigorously ir 
a future paper. 


We now summarize the search procedure for the system (2.6) where the 
N s -dimensional observation vector y" is given and an initial N q -dimensional 
parameter estimate 9 q is available. 

First the expected observation m and the covariance matrix R are determined from 
(3.3) and (3.S). The spectral decomposition of R as well as the quantities given by (b.17) 
- (3.33) then can be determined by standard "matrix algebra routines. And therefore 

from (3.32), (3.33) and (3.34) one obtains an N -dimensional gradient approximation g 

z 

and an N x N Her-ian approximation M. 

The parameter estimate 6 q can then be updated by making the quasi-Newton 
correction: 

6 t »9 - y M _1 g" (3.35) 

* o o ~ 

Here y q is an appropriate scalar chosen to improve the updated parameter 

estimate. In accordance w? i the geuera) theory of Ne'-rton iterations in function 
spaces [5], one can repeat this procedure until the solutions of the linearized problems 
converge to the solution of the underlying nonlinear problem. 

This analysis completes our outline of the maximum likelihood identification 
process. In Section 4 we give examples which illustrate the successful 
implementation of these schemes in useful applications. 
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•>. EXAMPLES: In this section we give examples of successful implementations of the 
previously discussed algorithms. 


We first consider the string (cf (2.13) 
Au = - (a(x) u'(x))' 


u(o) = u(l) = 0 


(4.1) 


For the case where the unknown tension parameter is constant, the point-source 
solution can be explicitly given: 

u(x|£) = * (1 - *>)(*<•> 

x^^ = max {x,f } (4.2) 

x < = min {x,f } 

And thus, as outlined in Section 2, all calculations could be given in terms of these 
quantities, without any truncation of the state space or the parameter space. 

In general, however, truncations in both spaces are necessary. For the string 
problem we consider an - dimensional state space of linear splines; the state 

variable then beccmes the vector of nodal values on the corresponding grid. For 
simplicity we take the grid to be uniform; thus, since the endpoints x = 0 and x = 1 
are fixed, we have: 


1 

Ax = 

N + 1 
x 

The state-space elements are then given by 
N 

r x (i) , , 

U = E U IC.(x) 
i=l 


where, as illustrated in Figure 3, the basis elements (r^x)} have the form 


! x -(i-l)Ax 
Ax 

x -(i)Ax 
Ax 


, (i-1) Ax < x < (i) Ax 


, (i) Ax < x < (i+1) Ax 


(4.3) 


(4.4) 


(4.5) 


0 


otherwise 




Pig. 3 Linear spline elements 

A similar discretization of the parameter space is possible. First we consider the 
augmented spline space 

N +1 

{r^xW W.6) 

1 UO 

where, as illustrated by Figure 3, the endpoint- elements r Q and + j are given by 

x 
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, 0 < x < Ax 


(Ax) - x 
Ax 

0 , A x < x < 1 

(4.7) 

0 , 0< x < (Ax) 

x - N x (Ax) 

, N (Ax) < x < 1 

x 

Ax 

Thus we have a corresponding parameter element 

N x +1 

a(x) = E a (l+1> r^x) (4.8) 

Uo 

which would give a parameter space {51 with dimension N^+2. 

However, as previously noted, the resolution of the parameter space often does 
not need to he as fine as the resolution of the state- space. We consider then the use 
of a piecewise linear parameter space of lower dimension where the only 
requirement is that the nodal points must be a subset of the nodal points of the 
state-space. The new parameter space is then a subset of the 2)- dimensional 

space given by (4.6). Let a be an Nq -dimensional parameter element (Nq S N^+2). 
Then a identifies with an element <T of the larger (N^ 2)- dimensional space and the 
relationship is given by 

a - £a (4.9) 

where B is an (N -t 2) x N Q matrix. And correspondingly, we have 

^ X V 

3a/3a = B . (4.10) 

This relationship simplifies the algorithms as described below since B is easy to 
construct, and the more cumbersome calculations which are needed to determine 
partial derivatives with respect to the parameter space are then specified in terms 
of the grid associated with the state space. Thus we have: 

a/aa-^B d/aa (4 . id 
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We illustrate these points with a sample calculation (see Fig. 4). We consider 
the case where there are seven sensors at the locations 

£ e {.125, .25, .375, .5, .625, .75, .35} (4.12) 


AAA 
S S S S S S S 

• 1 I 1 1 I I 1 . . « 

0.0 0.25 0.5 0.75 1 


S: SENSOR LOCATION 
A: ACTUATOR LOCATION 


rig. 4 String tension identification: sensor and actuator locations 
and three actuators at the locations 

£ e {.25, .5, .75}. (4.13) 

The data vector was derived from a plant with specifications 


plant 


(x) = 3 + x 


C f = 6(x- .25) +6(x-.5) + 6(x-.75) 


a = .001 
u 


(4.14) 


o = .001 
h 

For the state space we take the seven-dimensional space of linear splines {r^(x)} 

with nodes corresponding to the sensor locations (4.12), and for the parameter space 
we take the five -dimensional subset of linear splines with nodes corresponding to 
the set 


{0, .25, .5, .75, 1.} (4.15) 

The relaxation parameter y Q in (3.35) was chosen to speed up to the convergence of 

the iteration; these issues will be discussed more fully in a future report but we give 
the results of the calculations in Fig. 5. These numerical experiments appear to be 
very encouraging although with a crude approximation to the Hessian the 
convergence can be very slow. 
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ITERATE 

a* 11 

a‘ 2 » 

mm 

a< 4 > 

a ,5 > 

A 

INITIAL 

1.00 

a 

m 

1.00 

n 

B 

INTERMEDIATE 

3.19 

3.64 

1.98 

3.81 

3.62 

c 

FINAL 

2.98 

3.29 

3.55 

3.68 

4.05 

D 

TRUE 

3.00 

3.25 

3.50 

3.75 

Hi 


Fig. 5 Distributed parameter identification via spline analysis 


In a similar way, the algorithm was successfully applied to the wrap- rib antenna 
model (2.4). To simplify the calculations, we assumed here that the stiffness 
parameters were scalers although one could introduce a spline-based space as in the 
previous example. 

Again for simplicity we consider the case where there are six gores (N * 6), 
where a sensor is placed on the outer endpoint of each rib (r - 1), and where an 
actuator is placed at the midpoint of each rib (r *.5). This scheme is outlined in 
Figure 6. 


We introduce the set of N-dimensional unit vectors 



where the components of each are determined by 


(4.16) 
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• POINT - SENSOR LOCATION 

* POINT - ACTUATOR LOCATION 


Pig. 6 Rib stiffness identification: sensor and actuator locations 


1 , i = k 

0 , i*k 

The parameters of the likelihood functional are then given by 

a * .001 

a 

a = .001 
T1 

CT» £ «^6(r - .5) 

And the stiffness parameters of the pla^. are given by 
El = 1.25 .10* nt-m a 
T r = 1 ''S .10* 1 nt/m 
Tq 3 3.5 .10 nt/m 
L = 2.75 .10 -1 m 



(4.17) 

(4.18) 

(4.19a) 

(4.19b) 

(4.19c) 

(4.19d) 


We applied the algorithm then to the case where the unknown parameter was El 
while the other stiffness parameters were assumed to be known. 
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To discretize the suite- space eight equal subdivisions were made in the radial 
direction on each rib and in each mesh sector, in the circumferential direction five 
equal subdivisions were made in each mesh sector. The shape functions on the ribs 
were given by Hermite cubics while on the mesh the shape functions were given by 
splines linear in r and 9. In test cases this discretization produced at least three 
digits of accuracy in solving problems of the form (3.4). In all calculations the 
principle of cyclic symmetry (cf. (2.24)) was exploited to reduce the number of 
calculations. 

Convergence of the likelihood algorithm was very fast (see Pig. 7) when the 
relaxation parameter was taken to be 

y = 2.5 . 

'o 



0.100 0.125 0.160 0.175 

RIB STIFFNESS X 10* 6 ni-m 2 


ITERATION 
H\ -iBFR 
N 

ESTIMATED 

PARAMETER 

_ E, N 

RELATIVE ERROR. 

|< e i n -ei ES t)/ei est )| 

1 

1.75000* 10 6 

4. • 10* T 

2 

1.01930* 10 5 

2.-10' 1 

3 

1.20484* 10 5 

3. • 10^ 2 

4 

1.23866* 10 5 

2.*10" 4 

5 

1.23889* 10 5 

- 


ei plant “ 1 - 25000 * 1q5 
e, est ‘ 123889 * 1q5 


Pig. 7 Distributed parameter identification of beam stiffness parameter 
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Numerical experiments also demonstrated an improvement in the sensitivity of the 
identification schemes as the number of measurements was increased. Thus, Figure 
8 illustrates how, for the antenna problem considered, an increase in the number of 
sensors led to a steepening of the likelihood functional. Here the curves were 
shifted transversely for illustrative purposes. We note that no corresponding 
improvement in the parameter estimate occurred in these trials, possibly because of 
the less favorable signal- to- noise ratio which corresponds to sensing in the interior 
of the ribs. 


LIKELIHOOD FUNCTIONAL 



Fig. 8 Sensitivity of antenna stiffness - parameter identifica- 
tion according to number of aensors 

(a) One sensor per rib at r - 1.0 

(b) Two sensors per rib at r • 0.5, 1.0 

(c) Three sensors per rib at r - 0.5, 0.75, 1.0 

More detailed numerical experiments with distributed antenna stiffness 
parameters will be given in a future report. But the results outlined in this report 
demonstrate already the great potential for these algorithms. 
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